VLA Observations of DG Tau's Radio Jet: A highly coUimated thermal outflow 
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ABSTRACT 

The active young protostar DG Tau has an extended jet that has been weh studied 
at radio, optical, and X-ray wavelengths. We report sensitive new VLA full-polarization 
observations of the core and jet between 5 GHz and 8 GHz. Our high angular resolution 
observation at 8 GHz clearly shows an unpolarized inner jet with a size 42 AU (0.35^^) 
extending along a position angle similar to the optical-X ray outer jet. Using our 
nearly coeval 2012 VLA observations, we find a spectral index a — +0.46 =t 0.05, which 
combined with the lack of polarization, is consistent with bremsstrahlung (free-free) 
emission, with no evidence for a non-thermal coronal component. By identifying the 
end of the radio jet as the optical depth unity surface, and calculating the resulting 
emission measure, we find our radio results are in agreement with previous optical line 
studies of electron density and consequent mass-loss rate. We also detect a weak radio 
knot at 5 GHz located 7" from the base of the jet, coincident with the inner radio 



knot detected by Rodriguez et al. (2012) in 2009 but at lower surface brightness. We 
interpret this as due to expansion of post-shock ionized gas in the three years between 
observations. 



1. Introduction 



The evolution of a young stellar object (YSO) involves not only mass accretion, through a 
circumstellar molecular disk, but also the loss of angular momentum and mass flux through a 
narrowly colhmated jet, thought to be launched close to the YSO. Magnetic fields are suspected of 
playing an important role in both the launching and collimation of these jets ( Pudritz et al.|[20T2 



Cai et al. 2008). However, the processes through which these fields act are not well understood 



(e.g. Carrasco-Gonzalez et al.|[2010 ). 
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Classical T Tauri stars (CTTS) are YSOs which have evacuated most of their remnant cir- 
cumstellar envelope. This allows an unobscured line of sight to the central 100 AUs of their jet 
and measurements of forbidden line emission are possible ( Dougados et al.||2000|| ). Using the spec- 



troscopic diagnostic technique described by Bacciotti et al. (1999), the forbidden- line emission can 



provide constraints for the electron density, total hydrogen density, ionization fraction, and the aver- 
age excitation temperature of the gas in these jets. Consequently, probing the inner wind structure 
of CTTS is crucial in understanding the mechanisms by which these YSO jets are launched and 
collimated. 

Some CTTS have associated parsec-scale jets, with Herbig-Haro (HH) objects forming where 
the fast streams of material collide with the slower material along the jet ( [McGroarty et alT 2007 



McGroarty fc Rayl|2004 ). These jets are dynamical, evolving on the timescales of a few years, and 
contain knots which have proper motions in the range of a few 100 km s~^. In addition, about 
10 have been discovered to have X-ray emission. These X-rays most likely trace the fastest shocks 
in the jets of YSOs. If the material in the YSO jet is heated to X-ray emitting temperatures by 
shocks, shock velocities around 500 km s~^ required ( Schneider et al.pOll ). 



A large number of the compact jets associated with YSOs have been detected at radio wave- 
lengths. The dominant radio emission mechanism is thought to be thermal bremsstrahlung from 
the shock-heated gas. However, evidence for non-thermal emission from several YSOs has also been 
discovered. Polarization has been associated with several YSOs, including objects in the p Ophiuchi 
molecular cloud (Andre et al.'1988! White et al.|1992 ; Andre et al.|1992| ), the Taurus- Auriga molec- 
ular cloud ( [Phillips et aL][r993; Skinner ||1993trFeigelson et aLl|1994||Ray et al.||1997D , the R Coronae 
Australis region ( Choi et aL||2008 ) , the Orion region ( Zapata et al.||2004 ), and HH 7-11 (Rodriguez 
et al.||1999 ). Furthermore, a few YSOs have been associated with linearly polarized radio emission. 



these include HH 80-81 (Carrasco-Gonzalez et al. 2010), HD 283447 (Phillips et al. 1996), and 



the Orion Streamers ( Yusef-Zadeh et aL||1990 ). Additional characteristics of non-thermal emission, 
including strong variation in flux density on timescales of hours to days, a negative spectral index, 
and VLBI measurements of high brightness temperatures {Tb ^ 10^ K), have also been found 



(Curiel et al. 


1993; 


Hughes 1997; Andre et al. 


1992; Wilner et al. 19991 Rodriguez et al. 


2005) 


suggested source of this non-thermal emission is the gyrosynchrotron mechanism ( 


Andre||1996) 



Non-thermal emission is almost always associated with the more evolved weak-lined T Tauri 
stars (WTTS). This is consistent with the idea that the non-thermal coronal emission is revealed 
only after the optically thick mass outflows have evolved away dEisloffel et aLl|2000D . However, this 
idea may be too simple: Apparently non-thermal emission has recently been reported from several 
less-evolved YSO's that have infrared evidence for a disk ( Osten fc Wolk|[2009| ) . 



This paper focuses on DC Tau, a highly active CTTS driving a well studied energetic bipolar 
jet. Located in the Taurus Molecular Cloud (estimated distance 140 pc; Torres et al. ( 2009[ )), DC 
Tau is of spectral type K5-M0, with a mass of 0.67 M©, a luminosity of 1.7 Lq and an estimated age 
of 3 X 10^ yr ( [Kitamura et al.||l996b| |Gudel et all2007[ ). It was one of the first T Tauri stars to be 
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associated with an optical jet ( Mundt fc Fried|1983 ) and has been studied extensively with adaptive 
optics, interferometry and the Hubble Space Telescope. The "micro-jet", associated with the HH 



158 knot, extends out to 12" (e.g. Eisloffel &; Mundt 1998) at a position angle of 223° (Lavalley 



et al. 1997). Eisloffel fc Mundt] (1998) calculated the proper motion of 4 knots observed in the 
HH 158 jet located at distances between 2. 5^^-10. 0'^ ( 350-1400 au), velocities around 150 km s~^. 



Furthermore, Dougados et al. (2000) calculated the proper motions of knots located within 4.0'^ 
from the source of the HH 158 jet; they determine velocities of ^200 km s~^. 

The jet of DG Tau has an onion-like structure within 500 AU of the star, where faster, 
highly collimated gas is nested in wider slower material, with maximum bulk gas speeds reach- 
ing 500 km s~^ (Bacciotti et al. 2000; Lavalley- Fouquet et al. 2000). This velocity structure is 



expected if the jet material is launched from a range of disk radii ( Agra- Amboage et al. 2011). 



Moreover, for the slower jet material in the jet the velocity between the two sides of the jet is 
between 6-15 km s~^; this velocity shift may indicate that the jet is rotating ( Bacciotti et al.||2002 ). 
Using [Fe II] observations and averaging over the central 1.0^' , the mass-loss rate from the high 
velocity gas is determined to be (1.6 ± O.8)xlO~^M0 yr~^ and from the medium velocity gas 
(1.7 ± O.7)xlO~^M0 yr"-*^, giving a total mass-loss rate for the velocity range of 50 to 300 km 
of (3.3 =b l.l)xlO~^M0 yr~^. However, this value is a lower limit to the mass-loss rate from the 
atomic component of the jet since the [Fe II] emission does not probe the whole range of velocities 
seen at optical wavelengths (]Agra- Amboage et al. 2011). The mass-loss rate of DG Tau has also 



been estimated using [OI]A6300, giving a rate of 3xlO"'^M0 yr-^ ( [Hartigan et"aL]|l995| ). The most 



recent estimate for the mass loss through the jet atomic component is Maurri et al. (2012), with 
M = (1.2 ± 0.4) X 10~^ Mq yr~^. Near-infrared evidence for a counter-jet has been reported by 
Pyo et al.| ( [20031 ). The redshifted emission appears suddenly at -0.7^' which suggests that the inner 
part of the counter-jet is hidden behind an optically thick circumstellar disk. 

DG Tau shows strong millimeter continuum emission thought to arise from a compact dust 
disk around the star. Assuming a dust opacity coefficient /^^^^O. 02-0. 05 cm^ g~^ at 147 GHz, the 
spectrum is consistent with thermal emission from a disk having a radius of about 110 AU and a 



mass 0.01-0.06 Mq (Kitamura et al. 1996b). At larger scales, ^^CO observations have revealed a 



gas disk with radius of 2800 au, oriented with its major axis perpendicular to the jet of DG Tau 
( Kitamura et al.|[l996a ). 



The X-ray jet of DG Tau was first discovered by Giidel et al. (2005) using a Chandra ACIS-S 
observation, which showed very faint soft emission along the jet out to a distance of about 5.0'^ to 



the SW with position angle 225°. In addition, Giidel et al. (2005) found that DG Tau reveals a 



new type of X-ray spectrum that includes two emission components with vastly different absorbing 
column densities: A weakly attenuated soft spectral component associated with a plasma with 
electron temperatures of no more than a few MK and a strongly absorbed hard spectral component 
associated with a hot plasma of several tens of MK.The soft X-ray component is thought to originate 
at the base of the jet where the first shocks form, while the hard component is attributed to emission 



arising from a magnetospheric corona. Moreover, Schneider &; Schmitt (2008) demonstrated that 
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the soft and hard X-ray components have a separation of ^ 0.2^' and that the soft X-ray emission 
is coincident with emission from optical hnes indicating that this X-ray component is indeed from 
the jet. 

There are few previous radio observations of DG Tau. It was first detected in a VLA survey 



of the Taurus- Auriga region at 5 GHz (Cohen et al. 1982). Further observations at 15 GHz and 



1.5 GHz (Cohen &; Bieging 1986) found that the structure is elongated. Cohen &; Bieging (1986) 



suggested that the radio structure and spectrum is consistent with free-free emission from ionized 



gas in an outflowing jet. Rodriguez et al. (2012) reported radio knots in the extended jet at 



0.42^' and 6.98^' respectively. They determine that these radio components were coincident with 
previously detected optical knots. 

In this paper we present the results of a Jansky Very Large Array (VLA) multi-frequency 
campaign of DG Tau. The goals of these observations were to determine the morphology of the 
inner jet, search for possible non-thermal coronal radio emission, and investigate the nature of the 
radio knots at large distance from the star. 



2. VLA Observations 

The radio campaign comprised three observing epochs between June 2011 and April 2012. 
The June 18, 2011 observation used two 128 MHz bands centered on 8.33 GHz and 8.46 GHz in 
A configuration. In 2012 we conducted two C-array observations using the newly- available 2 GHz 
bandwidth capability of the VLA. The March 22, 2012 observation spanned the frequency range 
of 4.5 GHz to 6.5 GHz, while the April 15, 2012 observation spanned 7.9 GHz to 9.9 GHz. The 
details of these observations are listed in Table [H 

For all three observations both the receiver bandpass correction and the absolute flux density 



^The VLA is operated by the National Radio Astronomy Observatory, which is a facihty of the National Science 
Foundation operated under cooperative agreement by Associated Universities, Inc. 



Table 1. VLA Results 



Epoch 


Array 


Freq 


BW 


Time 


Total flux 


Peak flux 


Sky RMS 


Qfwhm 






(GHz) 


(GHz) 


(min) 


(mJy) 


(mJy/beam) 


(MJy) 


(arcsec) 


2000.83 


A 


43.3 


0.1 


240 


3.9 


0.6 


100 


0.05x0.04 


2011.46 


A 


8.5 


0.256 


198 


1.1 


0.4 


10.0 


0.38x0.26 


2012.22 


C 


5.4 


2.0 


29 


1.05 


1.02 


6.2 


4.3x4.0 


2012.29 


C 


8.5 


2.0 


30 


1.29 


1.16 


8.0 


2.7x2.4 
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scale was set using the amplitude calibrator 3C48. The phase calibration used the angularly nearby 
phase calibrator source J0403+2600. In addition, we made short observations of both 3C84 and 
3C138 in our June 2011 observation. These sources were used to solve for the antenna polarization 
leakage terms and linear polarization position angle. Data reduction and imaging was carried out 
using Common Astronomy Software Application (CASA) package distributed by the National Ra- 
dio Astronomy Observatory (NRAO). Additionally, we used the Astronomical Imaging Processing 
System (AIPS) to calibrate and map an archival 4-hour A-array VLA observation of DG Tau from 
November 1, 2000 at 43 GHz. 



3. Results 

3.1. Imaging the DG Tau Jet 

After flagging bad visibilities and applying the normal amplitude and phase calibration, clean 
maps in total intensity (Stokes /, see Figure [T]), circularly polarized intensity (Stokes V) and linear 
polarized intensity were produced using the standard CASA data calibration procedures. The peak 
and integrated flux density of both the total and polarized intensity was obtained using the CASA 
task IMFIT (see Table [T]). We did not detect either circular or linear polarization, with upper 
limits 0.03 mJy beam"-*^ (3%) for circular and 0.02 mJy beam"-*^ (2%) for linear polarization. Both 
the rising spectral index and lack of polarization suggest that any non-thermal coronal component 
is absent or very weak. 

The A-array 8 GHz total intensity contour map at epoch 2011.46 is shown in Figure [T] In 
this high-angular resolution image, the emission extends approximately 0.35^' SW in approximately 
same direction as the optical jet. Both the asymmetry and elongation of the emission suggest 



that the source is the collimated jet mapped by optical observations (Bacciotti et al. 2000; Agra- 



Amboage et al. 2011) which have resolved the inner jet to within 0.1^' of the central source and 



determined that collimation must occur within this region. 

The 2012 observations were taken in C-array, with much lower angular resolution but higher 
sensitivity to larger, low-surface brightness components. Contour maps of the 5 GHz (epoch 
2012.22) and 8 GHz (epoch 2012.29) observations are shown in Figure [2j The location of the 
central source for the 2012.22 (5.5 GHz) observation is coincident with the 2012.29 (8.5 GHz) cen- 
tral source within the centroid uncertainty (± 0.01'^). Note that the integrated flux density at 
8 GHz apparently increased 17% from epoch 2011.46 to 2012.29, but it is possible that the increase 
results from faint extended structure over-resolved by the smaller beam at epoch 2011.46 (cf. Table 
1). The 5.5 GHz map shows a component (labeled A) extending ^ 7^' SW whose centroid is coin- 
cident with a much weaker feature seen in the 8.5 GHz map. This component is nearly coincident 



with the knot component reported by Rodriguez et al. (2012) but at a much lower flux density; we 
discuss this in more detail in section 4.3. 
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Fig. 1. — : VLA 8 GHz naturally- weighted contour map of DG Tau's inner jet at epoch 2011.46. 
The contours are -3, 3, 6, 10, 15, 20, 25, 30, 35, and 40 x 10 /xJy beam~^, the RMS noise level of 
the image. The restoring beam is shown in the bottom left corner with dimensions given in Table [T] 

There is also a more distant feature in the 5 GHz map (labeled C) with a flux density S ^ 100 
/iJy about 14^' from the stellar position, but well-displaced from the jet axis. This feature is also 
weakly seen in the Rodriguez et al. map, and is certainly real. We do not know if it is associated 
with the DG Tau jet, although it is unlikely that an unrelated background source this strong would 
be located this close to the jet. We speculate that the feature may be a part of an extended bow 
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Fig. 2. — : VLA contour maps of DG Tau during the two spring 2012 observations. The 8.5 GHz 
map at epoch 2012.29 is shown on the right, while the 4.5-6.5 GHz emission (epoch 2012.22) map is 
on the left. Both maps were made used natural weighting with contour levels -3, 3, 4, 6, 8, 10, 12, 
15, 20, 40, and 60 x 8.0 and 6.2 /iJy beam~^, the RMS values of the respective maps. The restoring 
beam for each map is given in the bottom left corner with the dimensions given in Table [l} The 



cross indicates the location of knot A detected b> Rodriguez et al. (2012) in 2009 radio observations. 
Feature C, located ^20. 0^' SW of the inner centriod, is visible in both 2012 maps and on the 2009 



map of Rodriguez et al. (2012), but is displaced from the optical jet. 



shock associated with the optical knot seen near 13'^ ( Rodriguez et al.|[2012 ), but do not discuss it 
further in this paper. 



3.2. Spectral energy distribution 



Figure [3] shows a radio spectrum of DG Tau using our VLA data, previously published radio 
data (Cohen et al. 1982; Cohen fc Bieging| 1986; Rodriguez et al.| 2012| ), and an archival VLA 
observation of DG Tau at 43 GHz (project code AW545, PI David Wilner) that we cahbrated 
and mapped. The spectrum monotonically increases over the frequency range of 4.5 to 43 GHz. 
However, it is clear from the multiple observations at 8.5 GHz and 4.5 GHz that the flux density 
is significantly variable over a period of years. Therefore we use only our nearly coeval 2012 
observations to fit a power-law to the spectrum, giving a spectral index a— 0.46 =b 0.05. This 
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Fig. 3. — : Flux density of DG Tau as a function of frequency for observations at epochs 1981.22 
— 1985.34( [Cohen fc Bieging|[T986| ), 1994.29, 1996.98, and 2009.63 ( [Rodriguez eraLl[20T2l ) , 2000.83 
(archival VLA), and 2011.46 — 2012.29 (this paper). The solid line, a power-law fit to the nearly 
coeval 2012.22 and 2012.29 points only, has a spectral index a — 0.46 ± 0.05. 



spectral index is typical of collimated thermal jets, whose spectral indices are close to +0.6 for 
constant velocity isothermal outflows, but whose overall spectral index can vary from nearly flat to 
+2 depending on the physical conditions in the jet, such as velocity gradients and recombination 



in the flow ( [Reynolds 1986). To better constrain the spectral index of the radio emission we need 
nearly coeval observations over a much larger range of frequencies. 
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Fig. 4. — : Analytic functions used to model the electron density (left) and velocity (right) profiles 
across DG Tau's jet at 0.35^^ from the central source. The points are derived values of electron 
density and velocity from ( Maurri et al.||2012 ). 



4. Discussion 

4.1. Comparison between radio and optically-derived electron densities 

Free-free emission depends only on the plasma temperature, the observing frequency, and the 
linear emission measure i.e., the square of the electron density integrated along the line of sight 
to the observer. If the plasma temperature is known (e.g., from optical line observations), and 
the optical depth at a given location can be estimated (e.g., from source structure), the emission 
measure at that location can be estimated. The entire detected radio jet is the optically thick 
surface. We therefore identify the outermost detectable region of the radio jet with the unity 
optical depth surface (r = 1). We then can calculate the emission measure and compare with 
estimates of electron density and jet width at this location based on optical line ratios. 

Inspection of Figure 1 shows that the radio emission becomes undetectable ^ 0.35'^ (50 AU 
projected) from the base of the jet. Assuming a thermal jet, this should correspond to the r — 1 
surface, meaning inside ^ 0.35^' the jet is fully optically thick while outside this distance the jet 
is optically thin. The temperature of the jet gas is not well-constrained, but if the jet is launched 
via a quasi-steady centrifugal MHD disk wind, as suggested by the transverse velocity gradients 
( Bacciotti et al.||2002| ), the gas temperature should be not much higher than the photospheric 



temperature (4800 K). However, the innermost, high- velocity core is likely shock- heated and may 
have a temperature near 8000 K, typical to other stellar jets. The lower-speed, broader flow 
associated with warm molecular emission is cooler, with a temperature of 2000 K, determined from 



infrared line ratios (Takami et al. 2004). In the following we assume a mean temperature across 
the flow T = 5000 K. 
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The free-free optical depth can be written 
t{u,T,EM) = 1.06 ^ ^ ^ 



GH2 



104 K 



-1.35 



EM 



10^^ cm ^ 



(1) 



Assuming r = 1 at u— 8.5 GHz, the emission measure at 0.35" is 

EM = J Ue^ds = 3.3 X 10^*^ cm"^ 



(2) 



High angular resolution observations of optical line ratios have be used to determine both 



the electron density and jet width as a function of velocity bin along the jet (e.g., Coffey et al. 
2008; Maurri et al. |2012| ). By combining these measurements, we can estimate the axial electron 
density profile. Figure Qa) shows a simple Gaussian fit to the electron density as a function of axial 
distance, where the points represent the mean density in each velocity bin and the uncertainties are 



the FWHM width uncertainty in each bin, using data from Maurri et al. (2012). We have combined 



these measurements to fit a simple analytic function to the electron density as a function of jet 
axial distance p. 



where no = 2.5 x 10^ 
emission measure EM 



10^ cm ^ and 

l26 



ne{p) = no e v^-/ (3) 
5.5 au. Integrating this density profile, we find a linear 



3.2 X 10 cm ^, in excellent agreement with the radio data. 



4.2. Mass-loss rate of the ionized component 

We next use the axial electron density derived in the previous subsection to estimate the 
mass-loss rate of the ionized component of the jet outflow at r = 1 (50 AU projected distance). 
Figure |4]^b) shows a Gaussian fit to the FWHM widths of the velocity bins given in Maurri et al 



^2012D , 

V{p) = l/o e"(^) (4) 

where V^— 400 km s~^, and p^— 12 au. Using this velocity profile, along with the density profile 
determined in the previous sub-section, we determine the mass-flux of the ionized gas at 50 AU 
projected distance. 



poo 

/ V{p)n,{p) pdp. (5) 

^0 



M = 2 TT m7 

/o 

where rnl ^1.2 rrip is the average ion mass. Using the approximate analytic models shown in 
Figure [4] for ne(p) and V{p)^ we find mass flux M ^ 5 x 10~^ Mq yr~^, with more than half of 
the mass flux in the high- velocity component within 5 AU of the jet center. Since the high-velocity 
component is almost completely ionized ( Maurri et al.||2012 ), this mass-loss estimate should be less 



than a factor of two smaller than the mass-loss of both the ionized and neutral gas. 
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Recent mass loss estimates of the DG Tau jet include that of [Agra- Amboage et al. (2011) 



M = (3.3 ± 1.1) X 10-^ Mq yr-^ and 



Maurri et al. 



(2012), M = (1.2 ± 0.4) X 10-^ Mq yr' 



While these are both somewhat lower than the present model prediction, the differences are probably 
not significant given the poorly-constrained functional forms used for the axial density and velocity 
profiles. 



4.3. Radio Knots in the Extended Jet 



Radio observations of DG Tau at epoch 2009.6 (Rodriguez et al. 2012) show a radio knot 
located approximately 7'^ along the jet with an integrated flux density 150 /iJy, coincident with an 
optical [S II] knot. There is also a much weaker feature near 12^', also coincident with an optical 
knot. If these radio structures are associated with shock-compressed gas in the jet flow, we would 
expect them to evolve both in flux and position on the dynamical timescale of the shock and the 
expanding post-shock gas. 

Our observations confirm the existence of the inner knot. Figure [5] shows 5.4 GHz and 8.5 GHz 
images of DG Tau at epochs 2012.22 and 2012.29 respectively. At each frequency we show maps 
made with two different uv-plane density weighting functions: uniform weighting to maximize 
angular resolution (panels c, d), and natural (a.k.a. unity) weighting to maximize sensitivity to 
low- brightness features (panels a, b). The locations of the 7" knot (labeled A) and 12^' knot (labeled 
B) detected by Rodriguez et al.| ( |2012| ) are indicated by crosses. Both the 8.5 and 5.4 GHz natural 



weighted maps show radio emission coincident with knot A but not with knot B. We did not detect 
either knot A or B in our epoch 2011.46 A-array observation, probably because the knots were 
over-resolved, as discussed below. 



4.3.1. Knot proper motion 



The DG Tau jet optical knots move outward with sky-plane proper motions ^ 0.3''per year 
{V ^ 200 km/s), at least within 10 arcsec of the stellar position ( Eisloffel fc Mundt|1998 ; Dougados 



et al.| 2000| ). [Rodriguez et al. ( 2012| ) found radio knots that are approximately cospatial with 



the optical knots, and also move with this speed, suggesting that the radio and optical knots are 
two manifestations of the same traveling shocks. Here we examine whether our more recent knot 
observations are consistent with this hypothesis. We consider only the motion of knot A, since we 
did not detect knot B at any epoch. Note that all positions discussed below are angular separations 
from the stellar position in the sky plane. 
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Fig. 5.—: DG Tau contour maps at 8 GHz (a, c) and 5 GHz (b, d) at epochs 2012.29 and 2012.22 
respectively. The two maps in the first column of this figure, (a) & (b), are made using natural 
weighting. The contour levels are -3, 3, 4, 6, 10, 20, 40, 60, 80, 100, 120, and 140 times the RMS 
in each map; the RMS and beam dimensions are given in Table Q. The second column maps, 
(c) & (d), are made using uniform weighting in order to maximize angular resolution of the maps. 
The contour levels are -3, 3, 6, 10, 20, 40, and 60 times the RMS in each map. The RMS for (c) 
is 15 /iJy with beam dimensions 1.8'^ x 1.7^'; (d) has a RMS of 13 fiJj and beam dimensions of 
2.9^' X 2.6'^ The crosses in each of these maps indicate the locations of the two knots detected by 



Rodriguez et al. (2012). 
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Fig. 6. — : Projected knot location along jet vs. epoch (solid line) with observed centroid locations 
at epochs 2009.63 and 2012.22 
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Using the Rodriguez et al. (2012) calculated speed V = 198 km s ^ and their 2009.62 position 



for knot A, the expected position at epochs 2011.22 and 2011.29 is 7.8'^ We fit single-component 
GAUssian models to knot A in the 5.5 GHz image. At 5.5 GHz, knot A has a centroid position 
7.6'' zb 0.7'', while at 8.5 GHz it is 7.7" ± 0.3" (Figure [6]). Both positions agree with the predicted 
proper motion, but the fractional uncertainty is large, 37% in the higher resolution 8.5 GHz image. 

The optical position of knot A has recently been determined from HST STIS observations at 
epoch 2011.13 (Schneider et al. 2012, in preparation). Fitting a Gaussian to the spatial profile of 
the [S II] 673lA line within v= -265 =h 65 km/s results in a distance 7.2" it 0.1", compared with a 



predicted position 7.4" using the Rodriguez et al. (2012) proper motion. The measured position is 



only slightly offset from the predicted position, and could indicate that the peak radio and optical 
intensities occur along different lines of sight in the shocked emission region. This is plausible, 
since the radio and optical line emergent intensities have different dependencies on the density and 
temperature of the gas. 



4-3.2. Knot evolution: Expanding post- shock gas? 



In addition to moving with the predicted proper motion, the flux density of knot A decreased 
dramatically. In order to avoid differences associated with spectral index and angular resolution. 



we can compare the epoch 2009.63 measurement of Rodriguez et al. (2012) with our 2012.29 ob- 



servation, both of which were at 8.5 GHz and had the same angular resolution (VLA at C-array). 
The epoch 2009.63 integrated flux was (150 =b 20) /xJy while the 2012.29 integrated flux was 
(46 =b 26) /iJy, both measured using naturally-weighted maps. This is a 70% decrease in 2.7 years, 
and is certainly real given the matched frequency and angular resolutions of the two observations. 
This behavior is not unprecedented: Radio knots in several other YSO jets have been observed to 



decrease with time and become undetectable within years of their initial ejection (e.g. Marti et al 
1998]). The total flux density at 5.5 GHz (epoch 2012.22) was (73 ± 25) /xJy. This implies a 
spectral index —0.4 ± 0.4, consistent with optically thin thermal emission, but highly uncertain 
because of the large uncertainty in each measurement, and the differing angular resolutions. 



Rodriguez et al. (20121) proposed a model consisting of periodic generation of shocks moving 
outward at at projected speed 200 km s~^ in a conical outflow. The model predicts a periodic 
variation of knot flux density caused by corresponding periodic velocity variations in a bipolar 



outflow. Fig. 8 of Rodriguez et al. (2012) predicts a flux for knot A near that observed at epoch 



2012.29, but occurring 1.5 years after the flux level observed at epoch 2009.63 rather than 2.7 
years. However, since their model is parameterized by several geometrical parameters which are 
not well-constrained, a detailed comparison is not appropriate with only a single additional flux 
measurement. Instead, guided by Occam's razor, we interpret the flux decrease using a much simpler 
conceptual model of an expanding volume of post-shock gas characterized by a single variable, the 
electron density. 
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Fig. 7. — : (a) DG Tau knot A flux density vs. epoch calculated using optically thin expanding 
sphere model (see text), with observed values at epochs 2009.63 ( Rodriguez et al.|2012 ) and 2012.22 
(this paper), (b) Model electron density vs. projected axial distance (arc sec) at epochs 2009.63 
(solid line) and 2012.22 (dashed line), (c) Model surface brightness vs. epoch (solid line) with 
observed values (•) and upper limits (x). (d) Optical depth at 5 GHz (solid line) and 8 GHz 
(dashed line) vs. epoch, confirming the optically-thin model assumption. 
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We model the knot as an optically-thin, free-free emitting isotropic sphere whose density scales 
exponentially with radial distance p, 



neip) = 



exp 



-ln(2)( 

\PF 



(6) 



where 2pF is the fuh width at half maximum scale, and the normalization ensures that the integrated 
density is a constant i.e., the total number of electrons does not change as the sphere expands. 
We also assume the knot is isothermal, a simplistic but reasonable assumption since for optically- 
thin emission, the flux density dependence on temperature is much smaller than for density (S oc 
21-0.35^^2^ _ We assume T = 10^ K, based on optical hne ratio observations ( Agra-Amboage et al. 



2011). The sphere is assumed to expand linearly with time, p{t) = po + ~ ^o)- 



The flux and size of the sphere were calculated as a function of time and compared with 
observations. We solved for best-fit model parameters N = 4.2 x 10^^, po — 336 AU(2.4'' at d 
= 140 pc), and p = 21 AU yr~^ (0.15^^ Figure [t] shows the flux density, electron density 

profile, surface brightness, and optical depth as a function of epoch for the fitted model, along with 
observed values and upper limits. The agreement is well within the measurement uncertainty of 
the radio knot at all epochs, and has r <^ 1, as expected. 



Maurri et al. 



(2012) measure electron densities Ue ^ 10^ cm~^ along the jet at a projected 
distance near 4-5'^ For a conical flow with constant ionization fraction, the density scales as r~^, 
so we expect Ue ^ 500 cm~^ outside the radio knot at a projected distance of According to the 
model, the maximum knot density varied from 2000 to 1000 cm~^ from epochs 2009.6 - 2012.2, 
implying a density enhancement factor of 4 in 2009.6, decreasing to 2 in 2012.2. Assuming that 
the present expansion continues, the model predicts that the density contrast will vanish and that 
the knot will disappear within a few years. 



5. Summary 



We report multi-epoch VLA observations of the pre-main sequence star DG Tau's radio jet. 
The radio spectrum {a = 0.46 ± 0.05) and lack of polarization indicate that the emission is 
bremsstrahlung, with no evidence for a non-thermal coronal component. Assuming the end of the 
radio jet at 0.35^' (50 AU projected distance) is the r = 1 surface, we calculate the column emission 
measure. Assuming an 'onion skin' ionization model and azimuthal symmetry, we find that the 
mean electron density in the centre of the jet is n = 2.5x10^ cm~^. This agrees well with optical 



estimates at this location (Maurri et al. 2012). We model the electron density and velocity axial 
profiles to calculate the mass loss of the ionized component, Mion ^ 5 x 10~^ Mq yr~^, with 
more than half of the mass flux in the high-velocity component within 5 AU of the jet center. 
This mass loss is comparable to the total mass loss calculated using optical line observations (e.g., 
Agra-Amboage et al.pOll ; Maurri et al. 2012^, indicating that most of the mass loss in the jet at 
this location is in the ionized component. 
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We confirm the existence of a radio knot near 7" recently reported by Rodriguez et al. (2012). 
The knot proper motion is consistent with a projected speed V — 200 km s~^, as suggested by 



Rodriguez et al. (2012) and previous optical estimates. The flux density of the knot dramatically 



decreased between 2009.6 and 2012.2. We present a simple model for radio emission from the 
knot consisting of an optically-thin ionized sphere which expands linearly with time. The model 
predicts that the FWHM size of the radio knot increases from 340 AU (2.4") to 390 AU (2.75") 
from 2009.6 to 2012.2, while the central electron density decreases from 2000 to 1000 cm~"^. The 
resulting radio flux decreases from 150 ^Jy to 50 ^Jy, in agreement with observations. By scaling 
from previously published density measurements at closer distances along the jet, we find that the 
density enhancement factor of the knot decreased from 4 in 2009.6 to 2 in 2012.2, and that the 
knot will disappear completely within a few years. 
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